<!-- saved from url=(0014)about:internet -->
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8" />
<title>RTW Report - fy_.c</title>
<script language="JavaScript" type="text/javascript" src="rtwreport_utils.js"></script>
</head>
<body bgcolor="#eeeeee" text="#1122aa" onload="try {if (top) { if (top.rtwFileOnLoad) top.rtwFileOnLoad(window.document); else local_onload();}} catch(err) {};">
<p>
<table border="0" cellspacing="0" cellpadding="6" width="100%%" height="100%%"><tr><td width="100%%" valign="top" bgcolor="#ffffff">
<h4>File: <a href="../fy_.c" target="rtwreport_document_frame" id="linkToText_plain">fy_.c</a></h4>
<pre id="RTWcode">
<a name="1"><font color="#888888 "><i>    1   </i></font></a><font color="#117755"><i>/*</i></font>
<a name="2"><font color="#888888 "><i>    2   </i></font></a><font color="#117755"><I> * fy_.c</i></font>
<a name="3"><font color="#888888 "><i>    3   </i></font></a><font color="#117755"><I> *</i></font>
<a name="4"><font color="#888888 "><i>    4   </i></font></a><font color="#117755"><I> * Code generation for function 'fy_'</i></font>
<a name="5"><font color="#888888 "><i>    5   </i></font></a><font color="#117755"><I> *</i></font>
<a name="6"><font color="#888888 "><i>    6   </i></font></a><font color="#117755"><I> * C source code generated on: Mon Jun 22 23:07:02 2015</i></font>
<a name="7"><font color="#888888 "><i>    7   </i></font></a><font color="#117755"><I> *</i></font>
<a name="8"><font color="#888888 "><i>    8   </i></font></a><font color="#117755"><I> */</i></font>
<a name="9"><font color="#888888 "><i>    9   </i></font></a>
<a name="10"><font color="#888888 "><i>   10   </i></font></a><font color="#117755"><i>/* Include files */</i></font>
<a name="11"><font color="#888888 "><i>   11   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;rt_nonfinite.h&quot;</font>
<a name="12"><font color="#888888 "><i>   12   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;ode78_.h&quot;</font>
<a name="13"><font color="#888888 "><i>   13   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;fy_.h&quot;</font>
<a name="14"><font color="#888888 "><i>   14   </i></font></a><font color="#992211">#</font><font color="#992211"><b>include</b></font> <font color="#992211">&quot;ode78__data.h&quot;</font>
<a name="15"><font color="#888888 "><i>   15   </i></font></a>
<a name="16"><font color="#888888 "><i>   16   </i></font></a><font color="#117755"><i>/* Variable Definitions */</i></font>
<a name="17"><font color="#888888 "><i>   17   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_d_emlrtRSI">d_emlrtRSI</a> = <b>{</b> 41, <font color="#1122ff">&quot;fy_&quot;</font>,
<a name="18"><font color="#888888 "><i>   18   </i></font></a>  <font color="#1122ff">&quot;C:/\xe4\xba\x91\xe7\x9b\x98\xe5\x90\x8c\xe6\xad\xa5\xe6\x96\x87\xe4\xbb\xb6\xe5\xa4?/GTOC8/Lowthrust_mee_v0/fy_.m&quot;</font>
<a name="19"><font color="#888888 "><i>   19   </i></font></a><b>}</b>;
<a name="20"><font color="#888888 "><i>   20   </i></font></a>
<a name="21"><font color="#888888 "><i>   21   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_e_emlrtRSI">e_emlrtRSI</a> = <b>{</b> 52, <font color="#1122ff">&quot;fy_&quot;</font>,
<a name="22"><font color="#888888 "><i>   22   </i></font></a>  <font color="#1122ff">&quot;C:/\xe4\xba\x91\xe7\x9b\x98\xe5\x90\x8c\xe6\xad\xa5\xe6\x96\x87\xe4\xbb\xb6\xe5\xa4?/GTOC8/Lowthrust_mee_v0/fy_.m&quot;</font>
<a name="23"><font color="#888888 "><i>   23   </i></font></a><b>}</b>;
<a name="24"><font color="#888888 "><i>   24   </i></font></a>
<a name="25"><font color="#888888 "><i>   25   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_f_emlrtRSI">f_emlrtRSI</a> = <b>{</b> 85, <font color="#1122ff">&quot;fy_&quot;</font>,
<a name="26"><font color="#888888 "><i>   26   </i></font></a>  <font color="#1122ff">&quot;C:/\xe4\xba\x91\xe7\x9b\x98\xe5\x90\x8c\xe6\xad\xa5\xe6\x96\x87\xe4\xbb\xb6\xe5\xa4?/GTOC8/Lowthrust_mee_v0/fy_.m&quot;</font>
<a name="27"><font color="#888888 "><i>   27   </i></font></a><b>}</b>;
<a name="28"><font color="#888888 "><i>   28   </i></font></a>
<a name="29"><font color="#888888 "><i>   29   </i></font></a><font color="#112266"><b>static</b></font> emlrtRSInfo <a name="var_g_emlrtRSI">g_emlrtRSI</a> = <b>{</b> 14, <font color="#1122ff">&quot;sqrt&quot;</font>,
<a name="30"><font color="#888888 "><i>   30   </i></font></a>  <font color="#1122ff">&quot;C:/Program Files/MATLAB/R2013a/toolbox/eml/lib/matlab/elfun/sqrt.m&quot;</font> <b>}</b>;
<a name="31"><font color="#888888 "><i>   31   </i></font></a>
<a name="32"><font color="#888888 "><i>   32   </i></font></a><font color="#117755"><i>/* Function Declarations */</i></font>
<a name="33"><font color="#888888 "><i>   33   </i></font></a><font color="#112266"><b>static</b></font> <font color="#112266"><b>void</b></font> <a href="#fcn_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">eml_error</font></a>(<font color="#112266"><b>void</b></font>);
<a name="34"><font color="#888888 "><i>   34   </i></font></a>
<a name="35"><font color="#888888 "><i>   35   </i></font></a><font color="#117755"><i>/* Function Definitions */</i></font>
<a name="36"><font color="#888888 "><i>   36   </i></font></a><font color="#112266"><b>static</b></font> <font color="#112266"><b>void</b></font> <a name="fcn_eml_error">eml_error</a>(<font color="#112266"><b>void</b></font>)
<a name="37"><font color="#888888 "><i>   37   </i></font></a><b>{</b>
<a name="38"><font color="#888888 "><i>   38   </i></font></a>  <font color="#112266"><b>static</b></font> char_T cv0[4][1] = <b>{</b> <b>{</b> <font color="#1122ff">'s'</font> <b>}</b>, <b>{</b> <font color="#1122ff">'q'</font> <b>}</b>, <b>{</b> <font color="#1122ff">'r'</font> <b>}</b>, <b>{</b> <font color="#1122ff">'t'</font> <b>}</b> <b>}</b>;
<a name="39"><font color="#888888 "><i>   39   </i></font></a>
<a name="40"><font color="#888888 "><i>   40   </i></font></a>  emlrtPushRtStackR2012b(&amp;<a href="ode78__data_c.html#var_h_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_h_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">h_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="41"><font color="#888888 "><i>   41   </i></font></a>  emlrtErrorWithMessageIdR2012b(<a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>, &amp;<a href="ode78__data_c.html#var_emlrtRTEI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRTEI');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRTEI</font></a>,
<a name="42"><font color="#888888 "><i>   42   </i></font></a>    <font color="#1122ff">&quot;Coder:toolbox:ElFunDomainError&quot;</font>, 3, 4, 4, cv0);
<a name="43"><font color="#888888 "><i>   43   </i></font></a>  emlrtPopRtStackR2012b(&amp;<a href="ode78__data_c.html#var_h_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_h_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">h_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="44"><font color="#888888 "><i>   44   </i></font></a><b>}</b>
<a name="45"><font color="#888888 "><i>   45   </i></font></a>
<a name="46"><font color="#888888 "><i>   46   </i></font></a><font color="#112266"><b>void</b></font> <a name="fcn_fy_">fy_</a>(<font color="#112266"><b>const</b></font> real_T y[14], real_T auxdata_mu, real_T auxdata_Tmax, real_T
<a name="47"><font color="#888888 "><i>   47   </i></font></a>         auxdata_c, real_T ydot[14])
<a name="48"><font color="#888888 "><i>   48   </i></font></a><b>{</b>
<a name="49"><font color="#888888 "><i>   49   </i></font></a>  real_T W;
<a name="50"><font color="#888888 "><i>   50   </i></font></a>  real_T Z;
<a name="51"><font color="#888888 "><i>   51   </i></font></a>  real_T C;
<a name="52"><font color="#888888 "><i>   52   </i></font></a>  real_T b_y;
<a name="53"><font color="#888888 "><i>   53   </i></font></a>  real_T Ckep;
<a name="54"><font color="#888888 "><i>   54   </i></font></a>  real_T dv0[18];
<a name="55"><font color="#888888 "><i>   55   </i></font></a>  real_T c_y[6];
<a name="56"><font color="#888888 "><i>   56   </i></font></a>  int32_T i;
<a name="57"><font color="#888888 "><i>   57   </i></font></a>  int32_T i2;
<a name="58"><font color="#888888 "><i>   58   </i></font></a>  real_T B[18];
<a name="59"><font color="#888888 "><i>   59   </i></font></a>  real_T LtB[3];
<a name="60"><font color="#888888 "><i>   60   </i></font></a>  real_T LtB_mag;
<a name="61"><font color="#888888 "><i>   61   </i></font></a>  int32_T u;
<a name="62"><font color="#888888 "><i>   62   </i></font></a>  real_T ar;
<a name="63"><font color="#888888 "><i>   63   </i></font></a>  real_T au;
<a name="64"><font color="#888888 "><i>   64   </i></font></a>  real_T ah;
<a name="65"><font color="#888888 "><i>   65   </i></font></a>  real_T A_Bacc[6];
<a name="66"><font color="#888888 "><i>   66   </i></font></a>  real_T b_ar[3];
<a name="67"><font color="#888888 "><i>   67   </i></font></a>  real_T b_A_Bacc[6];
<a name="68"><font color="#888888 "><i>   68   </i></font></a>  real_T dW_dex;
<a name="69"><font color="#888888 "><i>   69   </i></font></a>  real_T dW_dey;
<a name="70"><font color="#888888 "><i>   70   </i></font></a>  real_T dW_dL;
<a name="71"><font color="#888888 "><i>   71   </i></font></a>  real_T dZ_dhx;
<a name="72"><font color="#888888 "><i>   72   </i></font></a>  real_T dZ_dhy;
<a name="73"><font color="#888888 "><i>   73   </i></font></a>  real_T dZ_dL;
<a name="74"><font color="#888888 "><i>   74   </i></font></a>  real_T dC_dhx;
<a name="75"><font color="#888888 "><i>   75   </i></font></a>  real_T dC_dhy;
<a name="76"><font color="#888888 "><i>   76   </i></font></a>  real_T d_y;
<a name="77"><font color="#888888 "><i>   77   </i></font></a>  real_T e_y;
<a name="78"><font color="#888888 "><i>   78   </i></font></a>  real_T f_y;
<a name="79"><font color="#888888 "><i>   79   </i></font></a>  real_T g_y[6];
<a name="80"><font color="#888888 "><i>   80   </i></font></a>  real_T dv1[36];
<a name="81"><font color="#888888 "><i>   81   </i></font></a>  real_T dv2[6];
<a name="82"><font color="#888888 "><i>   82   </i></font></a>  real_T b_Ckep[36];
<a name="83"><font color="#888888 "><i>   83   </i></font></a>
<a name="84"><font color="#888888 "><i>   84   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------ */</i></font>
<a name="85"><font color="#888888 "><i>   85   </i></font></a>  <font color="#117755"><i>/*  state dynamics */</i></font>
<a name="86"><font color="#888888 "><i>   86   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="87"><font color="#888888 "><i>   87   </i></font></a>  <font color="#117755"><i>/*  ydot = fy_(t, y , auxdata) */</i></font>
<a name="88"><font color="#888888 "><i>   88   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="89"><font color="#888888 "><i>   89   </i></font></a>  <font color="#117755"><i>/*  Input argumuents: */</i></font>
<a name="90"><font color="#888888 "><i>   90   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------------- */</i></font>
<a name="91"><font color="#888888 "><i>   91   </i></font></a>  <font color="#117755"><i>/*  t                    [1x1 double]          time                     [ TU]         */</i></font>
<a name="92"><font color="#888888 "><i>   92   </i></font></a>  <font color="#117755"><i>/*  y                   [14x1 double]          state variables     [  -  ] */</i></font>
<a name="93"><font color="#888888 "><i>   93   </i></font></a>  <font color="#117755"><i>/*  auxdata       [structure]               parameters         [  -  ] */</i></font>
<a name="94"><font color="#888888 "><i>   94   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="95"><font color="#888888 "><i>   95   </i></font></a>  <font color="#117755"><i>/*  Output argumuents: */</i></font>
<a name="96"><font color="#888888 "><i>   96   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------------- */</i></font>
<a name="97"><font color="#888888 "><i>   97   </i></font></a>  <font color="#117755"><i>/*  ydot              [7x1 double]          vector field           [  -  ] */</i></font>
<a name="98"><font color="#888888 "><i>   98   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="99"><font color="#888888 "><i>   99   </i></font></a>  <font color="#117755"><i>/*  External functions called: */</i></font>
<a name="100"><font color="#888888 "><i>  100   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------------- */</i></font>
<a name="101"><font color="#888888 "><i>  101   </i></font></a>  <font color="#117755"><i>/*  None */</i></font>
<a name="102"><font color="#888888 "><i>  102   </i></font></a>  <font color="#117755"><i>/*  */</i></font>
<a name="103"><font color="#888888 "><i>  103   </i></font></a>  <font color="#117755"><i>/*  Copyright (C) 2015/06/22 by Chen Zhang,  */</i></font>
<a name="104"><font color="#888888 "><i>  104   </i></font></a>  <font color="#117755"><i>/*  BeiHang University, chenzhang.buaa@gmail.com */</i></font>
<a name="105"><font color="#888888 "><i>  105   </i></font></a>  <font color="#117755"><i>/*  ------------------------------------------------------------------ */</i></font>
<a name="106"><font color="#888888 "><i>  106   </i></font></a>  <font color="#117755"><i>/*  -------------------------------------------------------- */</i></font>
<a name="107"><font color="#888888 "><i>  107   </i></font></a>  <font color="#117755"><i>/*  注意pp冲突 */</i></font>
<a name="108"><font color="#888888 "><i>  108   </i></font></a>  memset(&amp;ydot[0], 0, 14U * <font color="#112266"><b>sizeof</b></font>(real_T));
<a name="109"><font color="#888888 "><i>  109   </i></font></a>
<a name="110"><font color="#888888 "><i>  110   </i></font></a>  <font color="#117755"><i>/*  -------------------------------------------------------- */</i></font>
<a name="111"><font color="#888888 "><i>  111   </i></font></a>  W = (1.0 + y[1] * muDoubleScalarCos(y[5])) + y[2] * muDoubleScalarSin(y[5]);
<a name="112"><font color="#888888 "><i>  112   </i></font></a>  Z = y[3] * muDoubleScalarSin(y[5]) - y[4] * muDoubleScalarCos(y[5]);
<a name="113"><font color="#888888 "><i>  113   </i></font></a>  C = (1.0 + y[3] * y[3]) + y[4] * y[4];
<a name="114"><font color="#888888 "><i>  114   </i></font></a>  emlrtPushRtStackR2012b(&amp;<a href="#var_d_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_d_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">d_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="115"><font color="#888888 "><i>  115   </i></font></a>  b_y = y[0] / auxdata_mu;
<a name="116"><font color="#888888 "><i>  116   </i></font></a>  <font color="#112266"><b>if</b></font> (b_y &lt; 0.0) <b>{</b>
<a name="117"><font color="#888888 "><i>  117   </i></font></a>    emlrtPushRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="118"><font color="#888888 "><i>  118   </i></font></a>    <a href="#fcn_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">eml_error</font></a>();
<a name="119"><font color="#888888 "><i>  119   </i></font></a>    emlrtPopRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="120"><font color="#888888 "><i>  120   </i></font></a>  <b>}</b>
<a name="121"><font color="#888888 "><i>  121   </i></font></a>
<a name="122"><font color="#888888 "><i>  122   </i></font></a>  Ckep = muDoubleScalarSqrt(b_y) / W;
<a name="123"><font color="#888888 "><i>  123   </i></font></a>  emlrtPopRtStackR2012b(&amp;<a href="#var_d_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_d_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">d_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="124"><font color="#888888 "><i>  124   </i></font></a>  dv0[0] = 0.0;
<a name="125"><font color="#888888 "><i>  125   </i></font></a>  dv0[6] = 2.0 * y[0];
<a name="126"><font color="#888888 "><i>  126   </i></font></a>  dv0[12] = 0.0;
<a name="127"><font color="#888888 "><i>  127   </i></font></a>  dv0[1] = W * muDoubleScalarSin(y[5]);
<a name="128"><font color="#888888 "><i>  128   </i></font></a>  dv0[7] = (W + 1.0) * muDoubleScalarCos(y[5]) + y[1];
<a name="129"><font color="#888888 "><i>  129   </i></font></a>  dv0[13] = -Z * y[2];
<a name="130"><font color="#888888 "><i>  130   </i></font></a>  dv0[2] = -W * muDoubleScalarCos(y[5]);
<a name="131"><font color="#888888 "><i>  131   </i></font></a>  dv0[8] = (W + 1.0) * muDoubleScalarSin(y[5]) + y[2];
<a name="132"><font color="#888888 "><i>  132   </i></font></a>  dv0[14] = Z * y[1];
<a name="133"><font color="#888888 "><i>  133   </i></font></a>  dv0[3] = 0.0;
<a name="134"><font color="#888888 "><i>  134   </i></font></a>  dv0[9] = 0.0;
<a name="135"><font color="#888888 "><i>  135   </i></font></a>  dv0[15] = C * muDoubleScalarCos(y[5]) / 2.0;
<a name="136"><font color="#888888 "><i>  136   </i></font></a>  dv0[4] = 0.0;
<a name="137"><font color="#888888 "><i>  137   </i></font></a>  dv0[10] = 0.0;
<a name="138"><font color="#888888 "><i>  138   </i></font></a>  dv0[16] = C * muDoubleScalarSin(y[5]) / 2.0;
<a name="139"><font color="#888888 "><i>  139   </i></font></a>  dv0[5] = 0.0;
<a name="140"><font color="#888888 "><i>  140   </i></font></a>  dv0[11] = 0.0;
<a name="141"><font color="#888888 "><i>  141   </i></font></a>  dv0[17] = Z;
<a name="142"><font color="#888888 "><i>  142   </i></font></a>
<a name="143"><font color="#888888 "><i>  143   </i></font></a>  <font color="#117755"><i>/*  calculate | lambda' * B | */</i></font>
<a name="144"><font color="#888888 "><i>  144   </i></font></a>  c_y[0] = y[7];
<a name="145"><font color="#888888 "><i>  145   </i></font></a>  c_y[1] = y[8];
<a name="146"><font color="#888888 "><i>  146   </i></font></a>  c_y[2] = y[9];
<a name="147"><font color="#888888 "><i>  147   </i></font></a>  c_y[3] = y[10];
<a name="148"><font color="#888888 "><i>  148   </i></font></a>  c_y[4] = y[11];
<a name="149"><font color="#888888 "><i>  149   </i></font></a>  c_y[5] = y[12];
<a name="150"><font color="#888888 "><i>  150   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 3; i++) <b>{</b>
<a name="151"><font color="#888888 "><i>  151   </i></font></a>    <font color="#112266"><b>for</b></font> (i2 = 0; i2 &lt; 6; i2++) <b>{</b>
<a name="152"><font color="#888888 "><i>  152   </i></font></a>      B[i2 + 6 * i] = Ckep * dv0[i2 + 6 * i];
<a name="153"><font color="#888888 "><i>  153   </i></font></a>    <b>}</b>
<a name="154"><font color="#888888 "><i>  154   </i></font></a>
<a name="155"><font color="#888888 "><i>  155   </i></font></a>    LtB[i] = 0.0;
<a name="156"><font color="#888888 "><i>  156   </i></font></a>    <font color="#112266"><b>for</b></font> (i2 = 0; i2 &lt; 6; i2++) <b>{</b>
<a name="157"><font color="#888888 "><i>  157   </i></font></a>      LtB[i] += c_y[i2] * B[i2 + 6 * i];
<a name="158"><font color="#888888 "><i>  158   </i></font></a>    <b>}</b>
<a name="159"><font color="#888888 "><i>  159   </i></font></a>  <b>}</b>
<a name="160"><font color="#888888 "><i>  160   </i></font></a>
<a name="161"><font color="#888888 "><i>  161   </i></font></a>  emlrtPushRtStackR2012b(&amp;<a href="#var_e_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_e_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">e_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="162"><font color="#888888 "><i>  162   </i></font></a>  b_y = (LtB[0] * LtB[0] + LtB[1] * LtB[1]) + LtB[2] * LtB[2];
<a name="163"><font color="#888888 "><i>  163   </i></font></a>  <font color="#112266"><b>if</b></font> (b_y &lt; 0.0) <b>{</b>
<a name="164"><font color="#888888 "><i>  164   </i></font></a>    emlrtPushRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="165"><font color="#888888 "><i>  165   </i></font></a>    <a href="#fcn_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">eml_error</font></a>();
<a name="166"><font color="#888888 "><i>  166   </i></font></a>    emlrtPopRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="167"><font color="#888888 "><i>  167   </i></font></a>  <b>}</b>
<a name="168"><font color="#888888 "><i>  168   </i></font></a>
<a name="169"><font color="#888888 "><i>  169   </i></font></a>  LtB_mag = muDoubleScalarSqrt(b_y);
<a name="170"><font color="#888888 "><i>  170   </i></font></a>  emlrtPopRtStackR2012b(&amp;<a href="#var_e_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_e_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">e_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="171"><font color="#888888 "><i>  171   </i></font></a>
<a name="172"><font color="#888888 "><i>  172   </i></font></a>  <font color="#117755"><i>/*  optimal thrust direction by PMP */</i></font>
<a name="173"><font color="#888888 "><i>  173   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 3; i++) <b>{</b>
<a name="174"><font color="#888888 "><i>  174   </i></font></a>    LtB[i] = -LtB[i] / LtB_mag;
<a name="175"><font color="#888888 "><i>  175   </i></font></a>  <b>}</b>
<a name="176"><font color="#888888 "><i>  176   </i></font></a>
<a name="177"><font color="#888888 "><i>  177   </i></font></a>  <font color="#117755"><i>/*  optimal thrust magnitude by PMP */</i></font>
<a name="178"><font color="#888888 "><i>  178   </i></font></a>  <font color="#112266"><b>if</b></font> (-auxdata_Tmax * (LtB_mag / y[6] + y[13] / auxdata_c) &lt; 0.0) <b>{</b>
<a name="179"><font color="#888888 "><i>  179   </i></font></a>    <font color="#117755"><i>/*  Engine turn 'On' */</i></font>
<a name="180"><font color="#888888 "><i>  180   </i></font></a>    u = 1;
<a name="181"><font color="#888888 "><i>  181   </i></font></a>  <b>}</b> <font color="#112266"><b>else</b></font> <b>{</b>
<a name="182"><font color="#888888 "><i>  182   </i></font></a>    <font color="#117755"><i>/*  Engine turn 'Off' */</i></font>
<a name="183"><font color="#888888 "><i>  183   </i></font></a>    u = 0;
<a name="184"><font color="#888888 "><i>  184   </i></font></a>  <b>}</b>
<a name="185"><font color="#888888 "><i>  185   </i></font></a>
<a name="186"><font color="#888888 "><i>  186   </i></font></a>  <font color="#117755"><i>/*  optimal thrust acceleration */</i></font>
<a name="187"><font color="#888888 "><i>  187   </i></font></a>  ar = (real_T)u * auxdata_Tmax / y[6] * LtB[0];
<a name="188"><font color="#888888 "><i>  188   </i></font></a>  au = (real_T)u * auxdata_Tmax / y[6] * LtB[1];
<a name="189"><font color="#888888 "><i>  189   </i></font></a>  ah = (real_T)u * auxdata_Tmax / y[6] * LtB[2];
<a name="190"><font color="#888888 "><i>  190   </i></font></a>
<a name="191"><font color="#888888 "><i>  191   </i></font></a>  <font color="#117755"><i>/*  state dynamics */</i></font>
<a name="192"><font color="#888888 "><i>  192   </i></font></a>  A_Bacc[0] = 0.0;
<a name="193"><font color="#888888 "><i>  193   </i></font></a>  A_Bacc[1] = 0.0;
<a name="194"><font color="#888888 "><i>  194   </i></font></a>  A_Bacc[2] = 0.0;
<a name="195"><font color="#888888 "><i>  195   </i></font></a>  A_Bacc[3] = 0.0;
<a name="196"><font color="#888888 "><i>  196   </i></font></a>  A_Bacc[4] = 0.0;
<a name="197"><font color="#888888 "><i>  197   </i></font></a>  A_Bacc[5] = muDoubleScalarPower(W, 3.0) * auxdata_mu / (y[0] * y[0]);
<a name="198"><font color="#888888 "><i>  198   </i></font></a>  b_ar[0] = ar;
<a name="199"><font color="#888888 "><i>  199   </i></font></a>  b_ar[1] = au;
<a name="200"><font color="#888888 "><i>  200   </i></font></a>  b_ar[2] = ah;
<a name="201"><font color="#888888 "><i>  201   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 6; i++) <b>{</b>
<a name="202"><font color="#888888 "><i>  202   </i></font></a>    b_y = 0.0;
<a name="203"><font color="#888888 "><i>  203   </i></font></a>    <font color="#112266"><b>for</b></font> (i2 = 0; i2 &lt; 3; i2++) <b>{</b>
<a name="204"><font color="#888888 "><i>  204   </i></font></a>      b_y += B[i + 6 * i2] * b_ar[i2];
<a name="205"><font color="#888888 "><i>  205   </i></font></a>    <b>}</b>
<a name="206"><font color="#888888 "><i>  206   </i></font></a>
<a name="207"><font color="#888888 "><i>  207   </i></font></a>    b_A_Bacc[i] = Ckep * A_Bacc[i] + b_y;
<a name="208"><font color="#888888 "><i>  208   </i></font></a>  <b>}</b>
<a name="209"><font color="#888888 "><i>  209   </i></font></a>
<a name="210"><font color="#888888 "><i>  210   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 6; i++) <b>{</b>
<a name="211"><font color="#888888 "><i>  211   </i></font></a>    ydot[i] = b_A_Bacc[i];
<a name="212"><font color="#888888 "><i>  212   </i></font></a>  <b>}</b>
<a name="213"><font color="#888888 "><i>  213   </i></font></a>
<a name="214"><font color="#888888 "><i>  214   </i></font></a>  ydot[6] = -(real_T)u * auxdata_Tmax / auxdata_c;
<a name="215"><font color="#888888 "><i>  215   </i></font></a>
<a name="216"><font color="#888888 "><i>  216   </i></font></a>  <font color="#117755"><i>/*  --------------------------------------------------------- */</i></font>
<a name="217"><font color="#888888 "><i>  217   </i></font></a>  dW_dex = muDoubleScalarCos(y[5]);
<a name="218"><font color="#888888 "><i>  218   </i></font></a>  dW_dey = muDoubleScalarSin(y[5]);
<a name="219"><font color="#888888 "><i>  219   </i></font></a>  dW_dL = y[2] * muDoubleScalarCos(y[5]) - y[1] * muDoubleScalarSin(y[5]);
<a name="220"><font color="#888888 "><i>  220   </i></font></a>  dZ_dhx = muDoubleScalarSin(y[5]);
<a name="221"><font color="#888888 "><i>  221   </i></font></a>  dZ_dhy = -muDoubleScalarCos(y[5]);
<a name="222"><font color="#888888 "><i>  222   </i></font></a>  dZ_dL = y[3] * muDoubleScalarCos(y[5]) + y[4] * muDoubleScalarSin(y[5]);
<a name="223"><font color="#888888 "><i>  223   </i></font></a>  dC_dhx = 2.0 * y[3];
<a name="224"><font color="#888888 "><i>  224   </i></font></a>  dC_dhy = 2.0 * y[4];
<a name="225"><font color="#888888 "><i>  225   </i></font></a>  emlrtPushRtStackR2012b(&amp;<a href="#var_f_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_f_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">f_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="226"><font color="#888888 "><i>  226   </i></font></a>  b_y = y[0] * auxdata_mu;
<a name="227"><font color="#888888 "><i>  227   </i></font></a>  <font color="#112266"><b>if</b></font> (b_y &lt; 0.0) <b>{</b>
<a name="228"><font color="#888888 "><i>  228   </i></font></a>    emlrtPushRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="229"><font color="#888888 "><i>  229   </i></font></a>    <a href="#fcn_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">eml_error</font></a>();
<a name="230"><font color="#888888 "><i>  230   </i></font></a>    emlrtPopRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="231"><font color="#888888 "><i>  231   </i></font></a>  <b>}</b>
<a name="232"><font color="#888888 "><i>  232   </i></font></a>
<a name="233"><font color="#888888 "><i>  233   </i></font></a>  d_y = y[0] / auxdata_mu;
<a name="234"><font color="#888888 "><i>  234   </i></font></a>  <font color="#112266"><b>if</b></font> (d_y &lt; 0.0) <b>{</b>
<a name="235"><font color="#888888 "><i>  235   </i></font></a>    emlrtPushRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="236"><font color="#888888 "><i>  236   </i></font></a>    <a href="#fcn_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">eml_error</font></a>();
<a name="237"><font color="#888888 "><i>  237   </i></font></a>    emlrtPopRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="238"><font color="#888888 "><i>  238   </i></font></a>  <b>}</b>
<a name="239"><font color="#888888 "><i>  239   </i></font></a>
<a name="240"><font color="#888888 "><i>  240   </i></font></a>  e_y = y[0] / auxdata_mu;
<a name="241"><font color="#888888 "><i>  241   </i></font></a>  <font color="#112266"><b>if</b></font> (e_y &lt; 0.0) <b>{</b>
<a name="242"><font color="#888888 "><i>  242   </i></font></a>    emlrtPushRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="243"><font color="#888888 "><i>  243   </i></font></a>    <a href="#fcn_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">eml_error</font></a>();
<a name="244"><font color="#888888 "><i>  244   </i></font></a>    emlrtPopRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="245"><font color="#888888 "><i>  245   </i></font></a>  <b>}</b>
<a name="246"><font color="#888888 "><i>  246   </i></font></a>
<a name="247"><font color="#888888 "><i>  247   </i></font></a>  f_y = y[0] / auxdata_mu;
<a name="248"><font color="#888888 "><i>  248   </i></font></a>  <font color="#112266"><b>if</b></font> (f_y &lt; 0.0) <b>{</b>
<a name="249"><font color="#888888 "><i>  249   </i></font></a>    emlrtPushRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="250"><font color="#888888 "><i>  250   </i></font></a>    <a href="#fcn_eml_error" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'fcn_eml_error');" target="rtwreport_document_frame"><font color="#1122aa">eml_error</font></a>();
<a name="251"><font color="#888888 "><i>  251   </i></font></a>    emlrtPopRtStackR2012b(&amp;<a href="#var_g_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_g_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">g_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="252"><font color="#888888 "><i>  252   </i></font></a>  <b>}</b>
<a name="253"><font color="#888888 "><i>  253   </i></font></a>
<a name="254"><font color="#888888 "><i>  254   </i></font></a>  emlrtPopRtStackR2012b(&amp;<a href="#var_f_emlrtRSI" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_f_emlrtRSI');" target="rtwreport_document_frame"><font color="#1122aa">f_emlrtRSI</font></a>, <a href="ode78__mex_c.html#var_emlrtRootTLSGlobal" onclick="if (top) if (top.docHiliteMe) top.docHiliteMe(window, 'var_emlrtRootTLSGlobal');" target="rtwreport_document_frame"><font color="#1122aa">emlrtRootTLSGlobal</font></a>);
<a name="255"><font color="#888888 "><i>  255   </i></font></a>
<a name="256"><font color="#888888 "><i>  256   </i></font></a>  <font color="#117755"><i>/*  co-state dynamics */</i></font>
<a name="257"><font color="#888888 "><i>  257   </i></font></a>  g_y[0] = y[7];
<a name="258"><font color="#888888 "><i>  258   </i></font></a>  g_y[1] = y[8];
<a name="259"><font color="#888888 "><i>  259   </i></font></a>  g_y[2] = y[9];
<a name="260"><font color="#888888 "><i>  260   </i></font></a>  g_y[3] = y[10];
<a name="261"><font color="#888888 "><i>  261   </i></font></a>  g_y[4] = y[11];
<a name="262"><font color="#888888 "><i>  262   </i></font></a>  g_y[5] = y[12];
<a name="263"><font color="#888888 "><i>  263   </i></font></a>  dv1[0] = 2.0 * au;
<a name="264"><font color="#888888 "><i>  264   </i></font></a>  dv1[6] = 0.0;
<a name="265"><font color="#888888 "><i>  265   </i></font></a>  dv1[12] = 0.0;
<a name="266"><font color="#888888 "><i>  266   </i></font></a>  dv1[18] = 0.0;
<a name="267"><font color="#888888 "><i>  267   </i></font></a>  dv1[24] = 0.0;
<a name="268"><font color="#888888 "><i>  268   </i></font></a>  dv1[30] = 0.0;
<a name="269"><font color="#888888 "><i>  269   </i></font></a>  dv1[1] = 0.0;
<a name="270"><font color="#888888 "><i>  270   </i></font></a>  dv1[7] = dW_dex * muDoubleScalarSin(y[5]) * ar + (dW_dex * muDoubleScalarCos
<a name="271"><font color="#888888 "><i>  271   </i></font></a>    (y[5]) + 1.0) * au;
<a name="272"><font color="#888888 "><i>  272   </i></font></a>  dv1[13] = (dW_dey * muDoubleScalarSin(y[5]) * ar + dW_dey * muDoubleScalarCos
<a name="273"><font color="#888888 "><i>  273   </i></font></a>             (y[5]) * au) - Z * ah;
<a name="274"><font color="#888888 "><i>  274   </i></font></a>  dv1[19] = -dZ_dhx * y[2] * ah;
<a name="275"><font color="#888888 "><i>  275   </i></font></a>  dv1[25] = -dZ_dhy * y[2] * ah;
<a name="276"><font color="#888888 "><i>  276   </i></font></a>  dv1[31] = ((dW_dL * muDoubleScalarSin(y[5]) + W * muDoubleScalarCos(y[5])) *
<a name="277"><font color="#888888 "><i>  277   </i></font></a>             ar + ((dW_dL * muDoubleScalarCos(y[5]) - W * muDoubleScalarSin(y[5]))
<a name="278"><font color="#888888 "><i>  278   </i></font></a>                   - muDoubleScalarSin(y[5])) * au) - dZ_dL * y[2] * ah;
<a name="279"><font color="#888888 "><i>  279   </i></font></a>  dv1[2] = 0.0;
<a name="280"><font color="#888888 "><i>  280   </i></font></a>  dv1[8] = (-dW_dex * muDoubleScalarCos(y[5]) * ar + dW_dex * muDoubleScalarSin
<a name="281"><font color="#888888 "><i>  281   </i></font></a>            (y[5]) * au) + Z * ah;
<a name="282"><font color="#888888 "><i>  282   </i></font></a>  dv1[14] = -dW_dey * muDoubleScalarCos(y[5]) * ar + (dW_dey * muDoubleScalarSin
<a name="283"><font color="#888888 "><i>  283   </i></font></a>    (y[5]) + 1.0) * au;
<a name="284"><font color="#888888 "><i>  284   </i></font></a>  dv1[20] = dZ_dhx * y[1] * ah;
<a name="285"><font color="#888888 "><i>  285   </i></font></a>  dv1[26] = -muDoubleScalarCos(y[5]) * y[1] * ah;
<a name="286"><font color="#888888 "><i>  286   </i></font></a>  dv1[32] = ((-dW_dL * muDoubleScalarCos(y[5]) + W * muDoubleScalarSin(y[5])) *
<a name="287"><font color="#888888 "><i>  287   </i></font></a>             ar + ((dW_dL * muDoubleScalarSin(y[5]) + W * muDoubleScalarCos(y[5]))
<a name="288"><font color="#888888 "><i>  288   </i></font></a>                   + muDoubleScalarCos(y[5])) * au) + dZ_dL * y[1] * ah;
<a name="289"><font color="#888888 "><i>  289   </i></font></a>  dv1[3] = 0.0;
<a name="290"><font color="#888888 "><i>  290   </i></font></a>  dv1[9] = 0.0;
<a name="291"><font color="#888888 "><i>  291   </i></font></a>  dv1[15] = 0.0;
<a name="292"><font color="#888888 "><i>  292   </i></font></a>  dv1[21] = ah / 2.0 * muDoubleScalarCos(y[5]) * dC_dhx;
<a name="293"><font color="#888888 "><i>  293   </i></font></a>  dv1[27] = ah / 2.0 * muDoubleScalarCos(y[5]) * dC_dhy;
<a name="294"><font color="#888888 "><i>  294   </i></font></a>  dv1[33] = -ah / 2.0 * C * muDoubleScalarSin(y[5]);
<a name="295"><font color="#888888 "><i>  295   </i></font></a>  dv1[4] = 0.0;
<a name="296"><font color="#888888 "><i>  296   </i></font></a>  dv1[10] = 0.0;
<a name="297"><font color="#888888 "><i>  297   </i></font></a>  dv1[16] = 0.0;
<a name="298"><font color="#888888 "><i>  298   </i></font></a>  dv1[22] = ah / 2.0 * muDoubleScalarSin(y[5]) * dC_dhx;
<a name="299"><font color="#888888 "><i>  299   </i></font></a>  dv1[28] = ah / 2.0 * muDoubleScalarSin(y[5]) * dC_dhy;
<a name="300"><font color="#888888 "><i>  300   </i></font></a>  dv1[34] = ah / 2.0 * C * muDoubleScalarCos(y[5]);
<a name="301"><font color="#888888 "><i>  301   </i></font></a>  dv1[5] = -2.0 * auxdata_mu * muDoubleScalarPower(W, 3.0) / muDoubleScalarPower
<a name="302"><font color="#888888 "><i>  302   </i></font></a>    (y[0], 3.0);
<a name="303"><font color="#888888 "><i>  303   </i></font></a>  dv1[11] = auxdata_mu / (y[0] * y[0]) * 3.0 * (W * W) * dW_dex;
<a name="304"><font color="#888888 "><i>  304   </i></font></a>  dv1[17] = auxdata_mu / (y[0] * y[0]) * 3.0 * (W * W) * dW_dey;
<a name="305"><font color="#888888 "><i>  305   </i></font></a>  dv1[23] = dZ_dhx * ah;
<a name="306"><font color="#888888 "><i>  306   </i></font></a>  dv1[29] = dZ_dhy * ah;
<a name="307"><font color="#888888 "><i>  307   </i></font></a>  dv1[35] = auxdata_mu / (y[0] * y[0]) * 3.0 * (W * W) * dW_dL + dZ_dL * ah;
<a name="308"><font color="#888888 "><i>  308   </i></font></a>  dv2[0] = 1.0 / (2.0 * W * muDoubleScalarSqrt(b_y));
<a name="309"><font color="#888888 "><i>  309   </i></font></a>  dv2[1] = -muDoubleScalarSqrt(d_y) / (W * W) * dW_dex;
<a name="310"><font color="#888888 "><i>  310   </i></font></a>  dv2[2] = -muDoubleScalarSqrt(e_y) / (W * W) * dW_dey;
<a name="311"><font color="#888888 "><i>  311   </i></font></a>  dv2[3] = 0.0;
<a name="312"><font color="#888888 "><i>  312   </i></font></a>  dv2[4] = 0.0;
<a name="313"><font color="#888888 "><i>  313   </i></font></a>  dv2[5] = -muDoubleScalarSqrt(f_y) / (W * W) * dW_dL;
<a name="314"><font color="#888888 "><i>  314   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 6; i++) <b>{</b>
<a name="315"><font color="#888888 "><i>  315   </i></font></a>    A_Bacc[i] = b_A_Bacc[i] / Ckep;
<a name="316"><font color="#888888 "><i>  316   </i></font></a>    c_y[i] = -g_y[i];
<a name="317"><font color="#888888 "><i>  317   </i></font></a>  <b>}</b>
<a name="318"><font color="#888888 "><i>  318   </i></font></a>
<a name="319"><font color="#888888 "><i>  319   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 6; i++) <b>{</b>
<a name="320"><font color="#888888 "><i>  320   </i></font></a>    <font color="#112266"><b>for</b></font> (i2 = 0; i2 &lt; 6; i2++) <b>{</b>
<a name="321"><font color="#888888 "><i>  321   </i></font></a>      b_Ckep[i + 6 * i2] = Ckep * dv1[i + 6 * i2] + A_Bacc[i] * dv2[i2];
<a name="322"><font color="#888888 "><i>  322   </i></font></a>    <b>}</b>
<a name="323"><font color="#888888 "><i>  323   </i></font></a>  <b>}</b>
<a name="324"><font color="#888888 "><i>  324   </i></font></a>
<a name="325"><font color="#888888 "><i>  325   </i></font></a>  <font color="#112266"><b>for</b></font> (i = 0; i &lt; 6; i++) <b>{</b>
<a name="326"><font color="#888888 "><i>  326   </i></font></a>    g_y[i] = 0.0;
<a name="327"><font color="#888888 "><i>  327   </i></font></a>    <font color="#112266"><b>for</b></font> (i2 = 0; i2 &lt; 6; i2++) <b>{</b>
<a name="328"><font color="#888888 "><i>  328   </i></font></a>      g_y[i] += c_y[i2] * b_Ckep[i2 + 6 * i];
<a name="329"><font color="#888888 "><i>  329   </i></font></a>    <b>}</b>
<a name="330"><font color="#888888 "><i>  330   </i></font></a>
<a name="331"><font color="#888888 "><i>  331   </i></font></a>    ydot[7 + i] = g_y[i];
<a name="332"><font color="#888888 "><i>  332   </i></font></a>  <b>}</b>
<a name="333"><font color="#888888 "><i>  333   </i></font></a>
<a name="334"><font color="#888888 "><i>  334   </i></font></a>  ydot[13] = -(real_T)u * auxdata_Tmax * LtB_mag / (y[6] * y[6]);
<a name="335"><font color="#888888 "><i>  335   </i></font></a>
<a name="336"><font color="#888888 "><i>  336   </i></font></a>  <font color="#117755"><i>/* ---------------------------------------------------------------- */</i></font>
<a name="337"><font color="#888888 "><i>  337   </i></font></a><b>}</b>
<a name="338"><font color="#888888 "><i>  338   </i></font></a>
<a name="339"><font color="#888888 "><i>  339   </i></font></a><font color="#117755"><i>/* End of code generation (fy_.c) */</i></font>
<a name="340"><font color="#888888 "><i>  340   </i></font></a></pre>
</td></tr></table>
</p>
</body>
</html>